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ABSTRACT 

We report the discovery of two new large H i shells in the direction of the Coalsack nebula. Both 
shells were observed with the Parkes Radio Telescope as part of the Southern Galactic Plane Survey. 
The largest shell, GSH 304-00-12, is at a distance of ~ 1.2 kpc and has derived physical dimensions of 
280 X 200 pc. The second shell, GSH 305-1-01-24, is at a distance of ^ 2.2 kpc and has derived dimensions 
of 280 X 440 pc. We present a simple numerical model to show that GSH 305-1-01-24 most likely formed 
from stellar winds in the Centaurus OBI stellar association. There is associated radio, infrared and 
Ha continuum emission. Both shells are situated in the Sagittarius-Garina arm, with GSH 305-1-01-24 
more distant. The far edge of GSH 304-00-12 is at the near side of the arm and opens into the interarm 
region. We find no evidence for closure at the near side of the shell and therefore describe the geometry 
as conical. Emission from the near side of the shell may be lost in absorption by the Goalsack Nebula. 



Subject headings: ISM: structure, bubbles 
hues: ISM 

1. INTRODUCTION 



Galaxy: structure, kinematics and dynamics — radio 



H I shells and supershells have long been studied in ex- 
ternal galaxies, as well as in our own Milky Way, as a probe 
of the star formation feedback processes in the interstellar 
medium (ISM; e.g. Westerlund & Matthewson 1966; Heiles 
1979). These structures, which are often identified as voids 
in the interstellar neutral hydrogen (H i), range from tens 
to hundreds of parsecs in diameter. Most shells form from 
the combined effects of stellar winds and supernovae that 
inject large quantities of energy (10^"'^ — 10^^ ergs) into 
the ISM. A few studies have directly identified the shell 
power source as an OB association (e.g. Saken et al. 1992; 
Brown et al. 1995). In these cases, however, standard wind 
blown bubble models (e.g. Weaver et al. 1977) fail to prop- 
erly account for the shells' radii and expansion velocities. 
As noted by Oey & Massey (1995) the model typically 
predicts radii that are too large and expansion velocities 
that are too small when compared to the Large Magellanic 
Cloud (LMC) shell population. 

In the nearby Large and Small Magellanic Clouds the 
large-scale structure of the ISM is dominated by hundreds 
of expanding shells (Staveley-Smith et al. 1997; Kim et al. 
1998). Many of these have been studied in an attempt 
to understand the dynamics of shells (Oey 1996a, b; Points 
et al. 1999) and their role in the ISM. Unfortunately, much 
less is known about the dynamics, density and distribution 
of H I shells in the Milky Way. Early H i surveys of the 
Milky Way resulted in the discovery of many shells (Heiles 



1979, 1984). New high resolution surveys of the Galactic 
Plane, such as the Canadian and Southern Galactic Plane 
Surveys (Taylor 1999; McClurc-Grifhths et al. 2001) are 
increasing this number and our knowledge of the small- 
scale structure of Galactic H i. 

We here report on the identification of two large shells, 
GSH 304-00-12 and GSH 305+01-24, located at the edge 
of the nearby Sagittarius-Carina spiral arm. The Coalsack 
nebula lies between the Sun and the near edge of GSH 304- 
00-12, though it does not appear physically related to the 
shell. We believe that GSH 305-1-01-24 is associated with 
the Centaurus OBI stellar association. We employ a very 
simple model for a wind-blown bubble to show that the 
shell is well explained by the catalogued stars in Gen OBI. 
In §2 we present the observations and data analysis. In §3 
we describe the shells morphologically and dynamically 
and suggest a geometry to explain their characteristics. In 
§3.3 we explore the relationship of GSH 3054-01-24 with 
Cen OBI. 



2. OBSERVATIONS AND ANALYSIS 

The data presented here were obtained as part of the 
Southern Galactic Plane Survey (SGPS; McClure-Griffiths 
et al. 2001). These observations were made with the multi- 
beam system on the Parkes Radio Telescope, a 64 m an- 
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tenna situated near Parkes NSW, Australia^. The Parkes 
portion of the SGPS covers 253° < I < 358°, |6| = ±10°. 
The observations were made during four observing sessions 
on 1998 December 15-16, 1999 June 18-21, 1999 Septem- 
ber 18-27, and 2000 March 10-15. Observations were made 
by the process of mapping "on-the-fly" with the inner 
seven beams of the multibeam system. The telescope was 
scanned through three degrees in Galactic latitude while 
recording data in 5 s samples. 

These data wcirc; recorded in frequency switching mode 
using the narrow-band back-end system (Haynes et al. 
1999), with a total bandwidth at each frequency of 8 MHz 
across 2048 channels. The center frequency was switched 
between 1419 MHz and 1422.125 MHz every 5 s. To re- 
move the bandpass shape each sample was divided by the 
previous frequency switched sample and the residual band- 
pass shape fitted with a series of Fourier components. The 
spectra were then multipUed by the mean of the reference 
signal over the spectrum to reconstruct the emission with 
a flat baseline. This technique for bandpass calibration 
has the advantage that continuum information is retained 
but changes in the bandpass shape during the course of the 
observing day are removed. Absolute brightness temper- 
ature calibration of the H i line data was performed from 
daily observations of the lAU standard regions S6 and S9 
(Williams 1973). These data were shifted to the local stan- 
dard of rest (LSR) by applying a phase shift in the Fourier 
domain. A more detailed description of the Parkes observ- 
ing strategy and calibration is given by McClure-Griffiths 
et al. (2000). 

The calibrated data were imaged using Gridzilla, a grid- 
ding tool created for use with Parkes multibeam data and 
found in the ATNF package of A1PS++^. The gridding al- 
gorithm as implemented in Gridzilla is described by Barnes 
et al. (2001). The data presented here were gridded using 
a weighted median technique, assuming a beamwidth of 
16'. A Gaussian smoothing kernel of FWHM 18' was em- 
ployed with a cutoff radius of 10', and a cellsize of 4'. The 
gridded data were imported to the MIRIAD'' data reduc- 
tion package for analysis (Sault fc Killeen 2000). Off-line 
channels in the ranges +160 to +200 km s~^ and —145 to 

— 100 km s~^ were used for continuum subtraction with 
the MIRIAD task "IMLIN". The final, calibrated data 
have an angular resolution of 16', a velocity resolution of 
0.82 km s~^ and a rms noise of ~ 0.3 K. 

3. H I SHELLS GSH 304-00-12 AND GSH 305+01-24 

We here present images of a pair of new shells in the di- 
rection of the Southern Coalsack. The larger shell has an 
extremely large angular size of 29° x 20° and is centered 
at / = 303?9, b = -0?2, v = -11.5 km s'^. Following 
the lAU standard naming convention based on the cen- 
ter coordinates of the shell, we name this shell GSH 304- 
00-12. Velocity channel images of GSH 304-00-12 from 

— 14.8 km s~^ to —7.2 km s^^ in intervals of 2.5 km s~^ 
are shown in Fig. 1. The shell is a large elliptical structure 
extending to \b\ ~ 10° in an arc that is nearly symmetric 

' The Parkes telescope is part of the Australia Telescope which is funded by the Commonwealth of Australia for operation as a National Facility 
managed by CSIRO. 

® See http://aips2.nrax).edu/docs/aips++.html 

R. J. Sault and N. E. B. Killeen, 2000, The Miriad User's Guide (Sydney: Austraha Telescope National FaciUty), at 
http://www.atnf.csiro.au/computing/software/miriad/. 



about the plane. For clarity, the shell is shown again in 
Fig. 2, a composite images that includes an H i image of 
the shell at u = — 12 km s~^ (panel a), a latitude cut across 
the shell (panel b) , and a velocity profile through the shell 
(panel c). The H i image is marked with an ellipse that 
traces the edge of the shell. There is a vertical line marking 
the position of the latitude cut and a cross-hair marking 
the position of the velocity profile. 

The shell first appears near w « — 27 km s~^, embed- 
ded in the H I emission from the Sagittarius-Carina arm, 
where it is as an arc of length ~ 20° peaked at 6 « 9°. 
Between v = —15 km s~^ and v = —9 km s~^ the size of 
the shell is roughly constant. At the central velocity of 

V = —12 km s^^ the shell is brightest on the high longi- 
tude end of the images and traces a faint arc through the 
low-longitude end. Towards larger velocities (i.e. towards 
the local gas) it decreases in angular size and finally dis- 
appears at w « —4 km s"'^. The caps of the shell, on both 
the near and far sides are not readily detected. On the 
near side, at -y ~ —4 km emission from the shell is ob- 
scured by absorption in the Coalsack nebula and emission 
from the local gas. The geometry of this shell is best de- 
scribed as conical. If there is emission from the far cap at 

V > 30 km s~^, it is swamped by the H i emission from the 
Sagittarius-Carina spiral arm. As seen in Figs. 1 & 2, the 
emission void above the plane contains a vertical filament 
at I « 306°, extending from 6 « 3? 5 to the top of the shell. 
This filament appears to correspond to the Galactic worm 
GW 305.5+9.4 (Koo et al. 1992). 

Using the Galactic rotation model of Fich et al. (1989) 
and adopting the lAU standard values for the Sun's or- 
bital velocity, Oo = 220 km s^^, and Galactic center dis- 
tance, _Ro — 8.5 kpc, we estimate the distance to GSH 
304-00-12 from the central velocity oi v = —12 km s~^. 
Assuming a velocity error of ±10 km s~^ to account for 
streaming motions (Burton 1988), we determine a distance 
of £) = 1.2 ± 1.0 kpc. At this distance the shell has 
physical diameters of 280 pc and 200 pc in the longitude 
and latitude directions, respectively. The near edge of the 
Sagittarius-Carina arm is at 1.26 ±0.23 kpc (Seidensticker 
& Schmidt-Kaler 1989). Although the shell distance er- 
rors are large, the coincidence of the shell and spiral arm 
distances suggests that the shell is at the edge of the arm. 

GSH 304-00-12 has a velocity full width of Av = 
23 km s~^. Disentangling the spatial extent of the shell 
from the physical expansion is difficult. Because of Galac- 
tic rotation, structures that span hundreds of parsecs will 
be observed over a range of radial velocities, even if they 
are not expanding. The velocity-distance relationship can 
be partially separated by considering the velocity gradient 
along the line of sight. At the location of shell the velocity 
gradient is ~ 11 km kpc~^. Therefore, if the shell is 
roughly spherical, such that the shell diameter along the 
line of sight is on the order of 240 pc, it should extend 
over approximately 3 km s^^ in velocity space. Given 
the errors involved in deriving distances from the rota- 
tion curve, we believe that the velocity expanse due to the 
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shell's physical extent in negligible. Therefore, we esti- 
mate the shell's expansion velocity as half of the velocity 
full width, Uoxp ~ Au/2 = 12 kms~^. This expansion 
velocity is comparable to shells of similar size in the SMC 
and somewhat smaller than shells in the LMC (Kim 1998; 
Stanimirovic 1999). 

The second shell, GSH 305+01-24, has a smaller angu- 
lar size of 7° X 11?3 and is located at / = 305?1, b = +1°, 

V = —24 km s~^. This shell is centered slightly above the 
plane and expands further above and below the plane than 
it does along the plane. An H i image of GSH 305+01- 
24 at t; = —24 km is shown in Fig. 3 overlaid with 
the positions of stars in the Cen OBI stellar association. 
The relationship of the shell to the stellar association is 
discussed below. 

Like the first shell, GSH 305+01-24 does not have a 
clearly defined near cap. It does, however, seem to have a 
far cap. This shell appears to open conically from within 
the spiral arm at w ^ — 34 km s^^. It increases in an- 
gular size towards larger velocities in the interarm region 
and reaches its maximum size at t; = — 24 km s~^. The 
angular size is relatively constant until v ~ —20 km s~^, 
where the shell begins to close. Beyond f w 2 — km s"^, 
emission from the shell is confused with that of GSH 304- 
00-12. The central velocity of t; = —24 km implies a 
distance of £> = 2.2 ± 0.9 kpc. At this distance the shell 
has physical radii of 140 x 220 pc (Ixb). GSH 305+01-24 
has a velocity full width of Av = 14 km s^^. Again, the 
velocity width due to the physical extent of this shell along 
the line of sight is probably negligible so we calculate an 
expansion velocity of Vexp ~ 7 km s~^. 

3.1. System Geometry 

The ISM along this line of sight is complicated. Spectro- 
photometric mapping of the region show that the Coalsack 
nebula is at an average distance of 199 ± 40 pc and con- 
sists of at least two overlapping clouds at 188 and 243 pc 
(Seidensticker & Schmidt-Kaler 1989). CO emission from 
the Coalsack nebula is geometrically centered at I = 303°, 
6 = 0°, covers ~ 30 deg^ on the sky and extends at least 
over the systemic velocity range Vrad = —10 to +8 km s~^ 
(Bronfman et al. 1989). H i self-absorption from gas in 
the Coalsack was detected by Bowers et al. (1980) and 
is also observed in the SGPS dataset. The interarm re- 
gion between the Coalsack and the edge of the Sgr-Car 
arm is nearly devoid of dust and molecular gas and the 
near edge of the spiral arm is at about 1.3 kpc (Seiden- 
sticker & Schmidt-Kaler 1989). H i and CO emission be- 
tween I = 290° — 320° and covering LSR velocities between 
~ —40 and ^ —20 km s~^ corresponds to the Sgr-Car arm 
at a mean distance of ~ 2 kpc. For velocities in the range 

V « —30 km to w « —18 km s~^ the spiral arm is 
nearly perpendicular to the line of sight. There is a very 
slight velocity gradient, such that the high longitude end 
of the arm is at more negative velocities than the low lon- 
gitude end. 

For a small range of velocities (-27 <v< -20 km s"^) 
GSH 305+01-24 overlaps with GSH 304-00-12. Though 
the shells are coincident in velocity space, there is no rea- 
son to believe that they are physically overlapping. The 
central velocities place the shells more than a kiloparsec 

* See IRAS Sky Survey Atlas Explanatory Supplement at http://www. 



apart and the line of sight widths are insufScient to overlap 
the shells. Also, the morphology of the shells do not ap- 
pear altered by an interaction. The overlap is most likely 
due to the expansion of the shells and not a true interac- 
tion. 

Both shells appear to have formed in the edge of the 
Sagittarius-Carina arm. GSH 305+01-24 seems to close 
near the edge of the spiral arm, whereas GSH 304-00-12 ex- 
pands into the interarm region. We suggest that GSH 304- 
00-12 experienced exaggerated expansion into the interarm 
region due to the pressure gradient leading away edge of 
the spiral arm. This explains the relatively constant size 
over the range of velocities covering the interarm region 
and the lack of a detectable near cap. GSH 305-00-12 
also extends into the range of velocities covered by the 
Coalsack nebula. The nebula at a distance of 199 ± 40 
pc (Seidensticker & Schmidt-Kaler 1989), however, is con- 
siderably closer than GSH 305-00-12. Unless the distance 
to GSH 304-00-12 has been significantly over-estimated or 
the shell's diameter along the line of sight far exceeds its 
diameter in the plane of the sky, it is unlikely that the 
shell is physically interacting with the Coalsack. The lack 
of emission from a near cap, though, could alternately be 
explained by absorption by cold H i in the Coalsack. 

3.2. Related Structures 

The lower portion of GSH 305+01-24 was previously 
identified as an H i shell associated with the Wolf-Rayet 

star 9 Muscae (Cappa De Nicolau & Nicmela 1984). How- 
ever, given the larger field of view in the current observa- 
tions, this shell arc appears connected with the arc above 
the plane. The void of GSH 305+01-24 is bisected by emis- 
sion from the Galactic plane. Because the rotation curve 
is dual-valued at these velocities, some emission is likely 
due to much more distant gas located on the far side of 
the tangent point. In the context of these measurements 
the 9 Muscae shell must be part of the larger structure 
that we identify as GSH 305+01-24. Cappa De Nicolau 
& Niemela (1984) also noted a second H i depression at 
t; = -12 km s"^ that we identify as part of GSH 304-00- 
12. 

There is 2.4 GHz continuum emission corresponding to 
the lower half of GSH 305+01-24. Fig. 4 is a 2.4 GHz con- 
tinuum image from the Duncan et al. (1995) survey of the 
Southern Galactic Plane overlaid with contours of an H i 
channel image of the shell at u = —24 km s~^. The con- 
tinuum emission excellently traces the inside edge of the 
shell below the plane. Also, the continuum may trace a 
portion of the shell above the plane near I = 301°, though 
this emission is very weak making a definitive detection 
difficult. There is an arc of continuum emission through 
the the upper H i void that does not correspond to any 
H I features. These arcs were identified by Duncan et al. 
(1995) as a possible supernova remnant, G303.5+0. How- 
ever, close examination of IRAS^ 100 /im images of the 
region reveals a very faint infrared shell coincident with 
the radio continuum features. These features were also 
identified as an Ha shell, the Coalsack Loop, by Walker 
& Zealey (1998). The detection of the arcs in the infrared 
implies that the emission is thermal and therefore suggests 
that the arcs are not part of a supernova remnant. The 

ipac. caltech.edu / ipac /iras/issa. html 
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close morphological similarity of the lower arcs and the H i 
shell leaves little doubt that the two are associated. We 
point out that the Coalsack Loop is much larger in angu- 
lar size than the Coalsack nebula so that even though the 
Loop is behind the nebula, extinction in the nebula will 
not affect the Ha emission. We believe, therefore, that the 
continuum emission is the photo-ionized region interior to 
the H I shell. 

3.3. Is GSH 305+01-24 Associated with Centaurus OBI? 

The Cen OBI stellar association is centered at I = 304°, 
b = +0?5 with radial velocities in the range —38 km s^^ to 
— 10 km s~^. The adopted distance is 2.5 kpc (Humphreys 
1978), which puts the cluster, to within the errors, at the 
same distance as GSH 305-1-01-24. The positions of the 
stars are plotted on a channel image of GSH 305+01-24 at 
V = —22 km s^^ in Fig. 3. The morphology of the shell 
appears to trace the stellar distribution, suggesting an as- 
sociation. The shell extends further above the plane than 
below, a characteristic that is matched by the stellar distri- 
bution. The cluster contains 21 stars with measured lumi- 
nosities and temperatures and calculated ages and masses 
(Humphreys 1978; Kaltcheva & Georgiev 1994). The ear- 
liest stars in the association are of type 09, which suggests 
that more massive stars may have already exploded as su- 
pernovae. A Wolf-Rayet star, 9 Muscae, of spectral type 
WC6 is also in the cluster. This system is a binary, pos- 
sibly a triple, consisting of the W-R-l-0 close binary and 
a distant O supergiant (Hartkopf et al. 1999). A priori, it 
seems likely that this association is the energy source for 
the GSH 305+01-24. 

To explore whether GSH 305+01-24 could be explained 
as a wind- driven bubble associated with Cen OBI, we em- 
ployed a very simple one-dimensional numerical model to 
study the evolution of a shell powered by the stars in this 
association. The model follows standard equations for the 
evolution of a pressure-driven wind bubble (e.g. Weaver 
et al. 1977; Oey & Massey 1995). The system is simplified 
to contain a centrally located power source, in this case 
the OB association, that supplies a time-dependent wind 
luminosity which drives the expansion of a shell. We as- 
sume that the internal pressure and temperature are con- 
stant throughout the bubble interior. The system is fur- 
ther simplified by assuming an ideal, monatomic gas equa- 
tion of state, and a cool, homogeneous external medium. 
The bubble expansion can be described by three first-order 
cqiiations: the rate of expansion of the shell, the equa- 
tion of momentum conservation and the equation of en- 
ergy conservation. The rate of expansion, Vexp, is simply 
the first-derivative of the shell radius, Rgh, 

dRsh ^-^-^ 



dt 



= Ve 



and the momentum conservation equation, including the 
effects of external and internal pressure, Pg and Pi, respec- 
tively, is 

?) ' y'^) 



(P. 



dt Ms 

where Mg is the mass of the shell and pe is the density of 
the ambient medium. The equation of energy conservation 
is 

dEsh , „ dV 



where is the wind luminosity and V is the bubble vol- 
ume, V = (47r/3)i?^^. The internal energy of the shell, 
Egh is related to the internal pressure according to the 
relation 



Esh = 



7-1' 



(4) 



dt 



dt 



(3) 



where the adiabatic index is 7 — 5/3. The wind lumi- 
nosity, — is calculated from an estimate of 
the mass-loss rate, M, and escape velocity, Woo, of each 
star in the Cen OBI association. The external pressure 
is given by Pg = /3pC^7^^, where Cg is the local speed of 
sound in the external medium, which is assumed to be 
~ 10 km s~^ for a warm neutral medium at T ~ 5000 K. 
Using temperatures, luminosities, masses, and ages for 21 
Cen OBI stars tabulated by Kaltcheva & Georgiev (1994), 
we estimate the mass-loss rate from the empirically deter- 
mined mass-loss rates of OB stars made by de Jagcr et al. 
(1988). The terminal velocity is estimated from the OB 
stellar wind model of Leitherer et al. (1992). Using the 
ages of the stars we calculate the total wind luminosity 
as the sum of all stars that exist at each time step. Only 
stars for which masses, luminosities, and ages are avail- 
able are included in the model. The spectral type, escape 
velocities, mass-loss rates, and wind luminosities are given 
in Table 1. For simplicity, we assume a constant external 
density, — 1.67x 10^^^ gm cm~^, corresponding to one 
hydrogen atom per cubic centimeter. The model is calcu- 
lated with time steps of 1.5 x 10^ yr and the initial condi- 
tions are calculated from the self-similar solutions to the 
differential equations at the first time step, t = 1.5 x lO'' yr. 

The results of the numerical integration are shown in 
Fig. 5. We explored two cases: one where the luminosity 
increases with time as stars are "born" in the association 
and another case of constant luminosity equal to the sum 
of all stars, regardless of age. The time-dependent wind 
luminosity case is plotted with a solid line and the con- 
stant wind luminosity case is plotted with a dashed line. 
The vertical dotted line marks the age of the oldest star 
in the association (14.7 Myr), with time i = assumed 
to be the birth of that star. The expansion velocity fiat- 
tens near 14 Myr because of the high mass-loss rate of the 
Wolf-Rayet star, 6 Muscae. Several differences between 
the two models are of note. First, the shell radius for the 
constant wind luminosity case can be fit as Pgh cx t^/^ , the 
analytic self-similar solution to the problem. But, in the 
time-dependent wind luminosity case the radius goes as 
^3/5 initially and then is better fit as Rsh oc t. Second, the 
usual assumption of coeval star formation for stellar asso- 
ciations leads to a slight over-estimate of the total energy 
input in a wind blown bubble. This effect is most notice- 
able as an overall increase in final shell radius by a factor 
of ~ 1.4; expansion velocity is relatively insensitive to the 
total energy input for larger ages. It seems, then, that 
including the stellar ages in the wind luminosity model 
results a marginally better fit to the data. 

We find that GSH 305+01-24 is surprisingly well de- 
scribed by the Weaver et al. (1977) model of a wind- 
blown bubble with a time-dependent wind luminosity from 
Cen OBI stars. This model predicts a shell radius of 
^ 210 pc, which corresponds well with the measured radii 
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of 140 X 220 pc; whereas the time-independent luminos- 
ity model predicts a radius of ~ 300 pc. The model 
predicts an expansion velocity of ~ 11 km s~^ for both 
the time independent and dependent luminosity cases, 
which also agrees well with an inferred expansion veloc- 
ity of ~ 6 km s~^. The good agreement suggests that 
GSH 305+01-24 was indeed formed from the stellar winds 
of the Cen OBI association. The size and expansion ve- 
locity are consistent with a shell age approximately equal 
to the ago of the oldest star in the association 15 Myr). 
Kaltcheva & Georgiev (1994), however, derive an age of 
8.9 ± 3.1 Myr for Cen OBI, from the average of the popu- 
lation's ages. The stellar ages are evenly distributed from 
14.7 Myr to 3.6 Myr with no evidence for sequential star 
formation, but rather there is evidence for continuous star 
formation, which makes it difficult to define an age for the 
cluster. Using the wind luminosities and ages of the Cen 
OBI association stars wc find the total energy released in 
the form of stellar wind luminosity during the lifetime of 
the association is ~ 8.3 x 10^^ ergs. We can compare this 
to the expansion energy of the shell - the equivalent energy 
instantaneously deposited at the center of the shell - de- 
rived by Chevalier (1974) and given in equation 2 of Hcilcs 
(1979):' Ee = 5.3 x 1043„J.i2^i^4 ^3^i2_ account for the 
current radius and rate of expansion of GSH 305-1-01-24, 
the expansion energy is ~ 9.8 x 10®^ ergs, which agrees 
well with the total energy of the Cen OBI stellar winds. 

3.3.1. Limitations of the Model 

The model shows, as it was intended to do, that the 
origin of GSH 305-1-01-24 can be explained, with respect 
to the energetics, as stellar wind luminosity from the Cen 
OBI association. Despite the good reproduction of the 
shell characteristics by this model, it should be stressed 
that it is a over-simplification and cannot replace full mag- 
netohydrodynamical simulations of H I shells from OB as- 
sociations. No attempts have been made to include multi- 
dimensional modeling, radiative cooling, magnetic fields, 
or the spatial distribution of the energy input sources. Ad- 
ditionally, the model contains no morphological informa- 
tion that can be helpful when trying to fully understand 
shell dynamics. Here we summarize some of the compo- 
nents of a non-idealized ISM that have not been included 
in the calculations and how those might effect the results. 

The model neglects energy injection from supernovae 
because there is no direct evidence for any associated with 
this cluster. However, it is possible that there have been 
supernovae in the association. The effect of supernovae on 
the wind-driven bubble evolution is very similar to that of 
stellar winds as an input energy source (McCray & Kafatos 
1987; Chu & Mac Low 1990).' McCray & Kafatos (1987) 
show in their equation 1 that for a standard initial mass 
function, an association of A^* stars with mass greater than 
7 Mq , supernovae input a mean luminosity over the associ- 
ation lifetime of P « 6.3 x 10^''' ergs^^(iV^i?5i). However, 
only ~ 20% is available in the first 10 Myr. Therefore, 
if A^'^t ~ 15, supernovae would have increased the modeled 
radius and expansion velocity of GSH 305-1-01-24 by about 
a factor of two. 

As shown in Fig. 3 the stars of the Cen OBI association 
are distributed throughout the interior of GSH 305-(-01-24. 
The model, however, assumes all stars are at the center of 



the shell. The B stars may have migrated as much as 
~ 50 pc in the shell life-time (McCray & Kafatos 1987), 
which suggests that the stars could have been initially dis- 
tributed over a region < 100 pc in diameter. This inferred 
remaining distribution should be negligible on the calcu- 
lated shell radius and expansion velocity, but could affect 
the morphology of the shell. 

The shell morphology would likely be dominated by ex- 
ternal density variations that were not included in this 
model. As shown in studies of supershcU expansion in 
stratified media (e.g. Mac Low et al. 1989), shells ex- 
panding out of the plane of the Galaxy experience exag- 
gerated expansion and blow-outs in that direction. For 
GSH 305-1-01-24, the shell radius is on the order of the scale 
height of the H I disk, at which height the density gradient 
of the disk becomes significant. Also, as noted above, the 
shells are at the edge of the Sagittarius-Carina spiral arm, 
so the density and pressure should decrease into the inter- 
arm region. Both of these effects will undoubtedly lead to 
exaggerated expansion down both density gradients. Small 
scale density fluctuations and the development of instabili- 
ties in the course of shell expansion can also effect the shell 
shape and size and can lead to break-outs in localized re- 
gions. The effect of such break-outs would be to equalize 
the internal and external pressures and stall expansion. 

Finally, three dimensional simulations of shell expansion 
in multi-phased, magnetized, turbulent media show shells 
that arc far from circular and exhibit many merging cav- 
ities and tunnels rather than a fully empty void (Korpi 
et al. 1999). Uniform magnetic fields in the disk can act 
to temporarily confine shells to the plane, and thereby af- 
fect the shell radius (Tomisaka 1998). Clearly, the model 
employed here has numerous limitations, but it adequately 
shows that winds from the Cen OB 1 association can ac- 
count for the observed shell radius and expansion velocity. 

4. CONCLUSIONS 

We have shown images of two new H i shells, 
GSH 305-F01-24 and GSH 304-1-00-12, located at the edge 
of the Sagittarius-Carina spiral arm in the direction of the 
Coalsack nebula. Both shells are very large, with radii on 
the order of 200 pc. We employed a simple time dependent 
wind luminosity model based on the ages and luminosities 
of the stars in the Cen OBI association to describe the 
evolution of GSH 305+01-24. The model predictions agree 
well with the observed shell radius and expansion veloc- 
ity. We use this agreement, the positional coincidence of 
the shell and OB association, and the similarities between 
the shell morphology and the stellar distribution to argue 
that GSH 305+01-24 is a wind-blown bubble. This shell 
is associated with a 2.4 GHz radio, infrared and Ha shell 
that were previously posited to be an old supernova rem- 
nant (Duncan et al. 1995; Walker & Zealey 1998). Because 
of the infrared detection of the shell we believe that the 
emission is thermal and therefore not part of a supernova 
remnant, but instead emission from the shell itself. 

We also found a large angular diameter shell, GSH 304- 
00-12, that does not appear to have formed from any 
known stellar associations. This shell is embedded in the 
edge of the Sagittarius-Carina arm and expanded into the 
interarm region, resulting in a conical geometry. We find 
no corresponding emission features in other wavebands. 
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GSH 304-00-12 has no detectable near cap, but the ve- 
locities where the cap is expected overlap velocities in the 
Coalsack nebula. We believe that there is no physical asso- 
ciation between the nebula and the shell, but we suggest 
emission from the cap may be absorbed by cold H I in 
the nebula. Additionally, although GSH 305-1-01-24 and 
GSH 304-00-12 overlap in velocity space, there is no direct 
evidence for any interaction between them. 
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HI SHELLS BEHIND THE COALSACK 
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Fig. 1. — Four velocity images of the largest shell, GSH 304+00-12, behind the Coalsack. The images show velocity planes from —14.8 km 
to —7.2 km and are separated by 2.5 km s~^. The grey scale is linear from 10 - 90 K, as shown in the wedge at the right. The shell, 
centered at I = 303?9 and b = — 0?2, extends to 6 = ±10?5 at the center and is best seen as the arcs above and below the plane on the left of 
the images. 
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Fig. 2.— H l channel image at v = -12.1 km s"! of GSH 304+00-12 (a), latitude slice at / = 308?15 across the shell (b), and a velocity 
profile through the shell at I = 305? 43, b = 2?07 (c). The outline of the shell is marked with the black ellipse on the H i channel image. The 
vertical line marks the position of the latitude cut (panel b) and the cross marks the position of the velocity profile (panel c). 
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Fig. 3. — H I channel image at t> = —24 km s"-*^ of GSH 305+01-24 in tlie edge of the Sagittarius-Carina arm. The grey scale is linear from 
10 to 130 K, as shown in the wedge at the right. The shell, of angular diameter 7° x 11?3, is centered at I = 305? 1, b = -|-1?0. The stars mark 
the positions of stars in the Cen OBI association. (Humphreys 1978). 
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GLON (degrees) 

Fig. 4. — Grey scale image of the 2.4 GHz continuum emission of the GSH 305+01-24 region overlaid with H I contours at d = —24 km s~^. 
The grey scale is linear and goes from 0.1 to 0.9 mjy Bm~^ as shown in the wedge at the right. The continuum image is from the Duncan 
et al. (1995) survey of the Galactic Plane and has been filtered to enhance small scale structure as described therein. The H i contours are 
at 45, 50, 60, and 70 K of T;,. The faint arcs of continuum emission in the lower half of GSH 305+01-24 trace the interior edge H i shown as 
shown by the contours. 



HI SHELLS BEHIND THE COALSACK 11 



300 



200 



o 



100 




t (Myr) 



10 



15 



T 100 




Fig. 5. — Simple numerical model of the shell radius and velocity evolution calculated for a constant luminosity input equal to the sum of 
the Cen OBI stellar wind luminosities (dashed line) and for a time dependent total wind luminosity calculated from the ages, masses and 
luminosities of the Cen OBI population (solid line). The vertical dotted line marks the present age of the stellar association. The actual shell 
radius and expansion velocity are well modeled by a time-dependent wind luminosity. The flattening of the velocity curve near 14 Myr in the 
time dependent model is due to the Wolf-Rayet star. 



MCCLURE-GRIFFITHS ET AL 



Table 1 

Properties of massive stars in the Cen OBI stellar association included in the numerical model for GSH 

305+01-24. 



Star 


Spectral Type^ 


Ago 




Log M 










(Myr) 


(km s-^) 


(Mq yr-i) 


(10^= erg s"^) 


(10*8 erg) 


HD 110639 


Bl lb 


10.3 


1514.1 


-6.22 


4.35 


14.15 


HD 111613 


A2 lab 


11.9 


893.8 


-6.10 


2.00 


7.51 


HD 111904 


B9 la 


9.6 


979.5 


-5.30 


15.15 


45.91 


HD 111934 


B2 lb 


9.1 


1427.3 


-5.60 


16.13 


46.31 


HD 111973 


B5 la 


8.6 


1169.3 


-5.91 


5.30 


14.39 


HD 111990 


B3 lb 


11.0 


1341.4 


-5.40 


22.58 


78.37 


HD 112272 


BO. 5 la 


4.7 


1563.9 


-5.52 


23.28 


34.52 


HD 112366 


B9 la-Ia 


10.0 


1044.1 


-5.40 


13.68 


43.16 


HD 112364 


B0.5 lb 


12.0 


1607.1 


-5.90 


10.25 


38.81 


HD 112842 


B3 lb 


14.7 


1201.5 


-5.40 


18.11 


84.02 


HD 113012 


BO lb 


12.6 


1764.2 


-6.22 


5.91 


23.50 


HD 113422 


Bl la 


6.5 


1423.9 


-5.92 


7.68 


15.75 


HD 114011 


Bl lb 


12.9 


1637.7 


-6.14 


6.12 


24.92 


HD 114213 


Bl lb 


9.5 


1320.5 


-6.22 


3.31 


9.93 


HD 115704 


BO. 5 la-Iab 


3.7 


1568.0 


-6.15 


5.49 


6.40 


HD 116119 


AO la 


10.2 


1006.3 


-7.16 


0.22 


0.71 


HD 115363 


Bl la 


3.6 


1494.8 


-5.27 


37.82 


42.96 


HD 114122 


BO. 5 la-Iab 


7.6 


1560.5 


-6.15 


5.43 


13.03 


HD 114886 


09 II-III 


4.9 


2161.6 


-5.57 


39.64 


61.29 


HD 113421 


B0.5 III 


6.7 


2386.8 


-6.90 


2.26 


4.78 


HD 113708 


B0.5 III 


6.2 


2197.7 


-6.90 


1.92 


3.75 


WR 48"= 


WC6+09.5I 




2000.0 


-3.74 


2294.0 


217.17 



^Prom Humphreys (1978). 

^Total lifetime energy input, Etot = L^ty where t = age 

*^The W-R+09.5 close binary is assumed to exist in its current state of high mass loss for ~ 3 X 10^ yr. 
The mass loss rate and terminal velocity are from Leitherer et al. (1997). 



Note. — Only stars with observed temperatures, luminosities, masses and ages were used. Terminal 
velocities were estimated from the OB stellar wind model of Leitherer et al. (1992). The mass-loss rates 
are based on the empirically determined mass-loss rates of OB stars in de Jager et al. (1988). In the 
model, stars are "turned on" at 14.7 Myr minus their age. 



